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. FOREWORD 

The work described, in this report was performed at the Gr um m an Aerospace 
Corporation, Bethpage, New York, and administered by the Vibration Section of the 
Structures and Dynamics Division, NASA Langley Research Center, Hampton, 
Virginia, 

The work performed under NASA Contract Nx\Sl-10635-21 with supplementary 
funding provided by the Space Division, Rockwell International (POM3WXMZ-4S3002) 
included development of a fondamental finite element hydroelastic formulation appli- 
cable to NASTRAN, implementation of the theoretical developments into NASTRAN 
and verification, and demonstration of the new technique on various problems in- 
cluding the 1/8-scale space shuttle external tanlc model. 
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ABSTRACT 

A fundamental reformulation of the NASTRAN hydroelastic analysis is 
developed in Volume I on the basis of Toupin's complementary variational principle 
of classical mechanics and on the basis of a physical interpretation of the NASTRAN 
fluid matrix equations. Emphasis is placed on the special case of an incompressible 
fluid model which is applicable to propellant tank hydroelastic analysis. A concise 
fluid inertia representation results from the assumption of incompressibility and 
the NASTRAN hydroelastic equations reduce to a simplified form associated with 
non-fluid structures, The efficiency of the incompressible hydroelastic formulation 
is enhanced for both fluid and structure by introduction of harmonic reduction as an 
alternative to Guyan redviction. The theoretical developments are implemented in 
NASTRAN by use of ALTER-DMAP statements and the modified NASTRAN hydro- 
elastic analysis technique is verified and demonstrated as an efficient and accurate 
approach by several illustrative problems. 

In Volume li the detailed modifications to the NASTRAN hydroelastic analysis 
are described. The DMAP operations are listed and applied to the l/8-Scale 
Shuttle Model External Tank. The resulting mode shapes and frequencies are 
compared with experimental results. 
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1 - INTRODUCTION 

A numerieally efficleat finite element hydroelastic analysis technique and a 
special, purpose redui^tion technique for geometrically axisymmetric structures have 
been developed and preseuti^d in Volume I (low number NASA CR to be published). 

While the formulation is quite general, it is particularly applicable to the NASTRAN 
program which contains a set of fluid finite elements currentl^r limited to axisymmetric 
geometry. The present report contains details of the implementation of the new for- 
mulation into NASTRAN including descriptions of the DMAP statements required for 
conversion of the program and details pertaining to problem definition and bulk data 
considerations. Details of the current l/S-scale space shuttle external tank mathe- 
matical model, numerical results and analysis/test comparisons are also presented. 
The appendices include a description and listing of a FORTRAN program used to 
develop harmonic transformation bulk data (multipoint constraint statements) and 
sample bulk data information for a number of hydroelastic problems. 
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2 - MODIFICATION OF THE NASTEAN HYDROELASTIC ANALYSIS 

Modification of the NASTEAN hydroelastic analysis according to the theoretical 
considerations discussed -Dreviously has been implemented by the use of ALTER- 
DMAP statements in Rigid Formats 7 {complex eigenvalue analysis), 3 (normal modes 
analysis), and 13 (undamped vibration analysis "with static pre-load effects).* Rigid 
Format 7 has been modified and condensed considerably since its function in the 
present application is to calculate fluid matrix data. The fluid matrix data formed in 
Rigid Format 7 is output onto cards or tape and used as input data for a normal modes 
analysis in either Rigid Format 3 or 13. Rigid Format 3 is used in cases where 
differential stiffness due to pressurization and/or fluid weight is negligilile and Rigid 
Format 13 is used when such effects are significant. The modifications in Rigid 
Formats 3 and 13 are minimal, consisting of the addition of the fluid mass matrix to 
the structural mass matrix in the anal^^sis set and the recovery of pressure and 
dependent displacement data after calculation of the modal solutions. In addition to 
the modifications of the Rigid Formats, special bulk data information must be formed 
to perform the operations required in the present applications. 

2. 1 DESCRIPTION OF DMAP OPERATIONS FOR FORMULATION OF FLUID DATA 

The following is an outline of the major operations performed in Rigid Format 7 
for the preparation of fluid matrix data. Operations utilized in both the general case 
(gravitational and/or ullage stiffness included) and the special case (zero free surface 
pressure) are presented. 

3. GPl generates coordinate system transformation matrices, tables of 
grid point locations, and tables for relating internal and external grid 
point numbers. 

8. GP2 generates Element Connection Table w’ith Internal indices. 

23. TAl generates element tables for use in matrix assembly and stress 
recovery. 


*Note: Rigid Format 13, although not documented in the latest NASTRAN 
manuals is included in the level 15 series programs. 
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28. SMAl generates stiffness matrix [k^ 1 , structural damping matrix 

r 4 1 ^ 

L ggJ Point singularity Table. 

42. Equivalence [i^g] to no general elements . 

44. Go to DMA P No. 47 if no general elements. 

45. SMA J adds general elements to [i*^g] to obtain stiffness matrix [i^gg] • 

49. GP4 generates flags defining members of various displacement sets 
(USET) and forms multipoint constraint equations 

52. Equivalence [k ] to [k ] , [m ] to [m ] , 


r 4 1 r 4 

to if no multipoint constraints. 

54. Go to DMAP No. 59 if general elements present. 

55. Go to DMAP No. 5 9 if no structural elements. 

5G. GPSP determines if possible grid point singularities remain. 

57. OFP formats table of possible grid point singularities and places it on the 
system output file for printing. 

60. Go to DMAP No. 65 if MCEl and MCE2 have already been executed for 
current set of multipoint constraints. 

61. MCEl partitions multipoint constraint equations [ftg] = [^ml^nl 

solves for multipoint constraint transformation matrix J 

N- 

63. MC E2 partitions stiffness, matrix 


gg" 


nn-^ 


R 


;] K} = 


- 0 . 


[b ] to [b 1 and 
gg-* nn-' 


1-1 


K = 
gg 


K K 
nn nm 


K K 
mn mm 


66 . 


and performs matrix reduction 

[k 1 - [k 1 [gH [k ] [k^ ] [g ] [g^] [k 1 [g 1, 

t nnJ t nnJ L mJ L mni L mnJ L mJ L m J L mmJ L mJ ’ 

Equivalence [k^] , if no single-point constraints. 
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68. Go to DMA P No. 71 if no single-point constraints. 

69. SCEl partitions out single-point constraints 


K = 
nn 


K, 


ff 


Kx- 

ts 


K r K 
sf ss 


72. Equivalence K... to K if no omitted coordinates 

^ IT ^ 3.3. 


74. Go to D1VL4P No. 88 if no omitted coordinates. 

95. BMG generates DMIG card images describing the interconnection of the 
fluid and the structure. 


100. MTRXIN generates fluid boundary matrices Aj^ and if a 

fluid structure interface is defined. The matrix is generated 

only for a nonzero gravity in the fluid. 

104. Go to next DhlAP instruction if cold start or modified restart. LBL13 
will be altered by the Executive System to the proper location inside the 
loop for unmodified starts within the loop. 

105. Beginning of loop for additional sets of direct input matrices. 

107. CASE extracts user re from CASECC for current loop. 


110. MTRXIN selects the direct input matrices for current loop, 

L irir J 

114. Equivalence [^^pp] [^*^ppj fi] 

115 . GotoDMAPNo. 112 if no [k,^ fj)J - 

116. Equivalence to [l^pp] if L^^pp] ' 

117. Go to DMAP No. 112 if no J ] . 

I, J 

118. ADD assembles matrix [f^p] ~ f^] [^"^pp] ' 

122. ADD subtracts [Aj^ ] from [Kp^j to obtain . 
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130. Equivalence to if no constraints applied. 

The following set of statements is utilized when either gravitational stiffness is 
ineludedi in the analysis or free surface displacement data is desired. Copies of the 
ALTER-DMAP statements and the complete DMiVP listing are presented Ln Figs. 2-1 
and 2-2, respectively 


136 a. Partition 

area matrix 


K1 


to form the fluid inertance matrix 


“ K] ’ 


the surface 

[■^em] surface gravitational stiffness matrix 

% ^ 

A, 


[^dd> 


em 


K 


The vector partitions the "d" set as follows 

{ 1 } 


136h. Partition to form free surface and structural interface matrices 

(Note: in the case where gravitational effects are included the free 
surface pressure subset consists of a single zeroth harmonic pressure) 




136c. Partition [a 1 to form free surface and structural interface matrices 
L emj 

(Note: in the case where gravitational effects are included the free 
surface displacement subset consists of a single zeroth harmonic 
displacement) 


'Si 

'^12 

Si 

'^2 


A 

em 


A 


Im 


‘2m 




.J 


y 


2-4 


id 


.... ii 

K t 

•wwvk-i 



Volume H 


136d Partition K to form free surface and structural interface matrices 
e 



136e. Calculate the transpose of 

136f. Solve for the pressure deviation recovery matrix 

hw] “ ^ 2 ]'" [^L] 

136g. Calculate the symmetric fluid mass matrix 

[^fLd] "" [^2m] Pdu2] 

136k, 1. solve [m^]^ [ 1 ^ 12 ] 

and form the displacement recovery matrix 

[’^ul2] " [^ 2 ] du2] 

136m. Output onto cards or tape: > [^du2] ’ ['^ul2] * 

136n. Form the reduced gravity stiffness matrix 

W= [w]Kii] kJ " K22] 

136o. Solve for the pressure recovery matrix (for the case of nonzero gravita- 
tional stiffness) 

[vl'-h2^0W 

136p. Output into cards or tape 

[■^eb]- kl]- M 

170. EXIT 

*OXJTPUT 3 for punched card PMC output is utilized in the listings presented in 
this report; if tape output is desired OUTPUT 1 should be used. 

2-5 



Volume H 


N 

S 


2 * 
* Z 

U s 


Ul 

> 


IL It 
U o 

> V 
V X 

IL It 
U. IL 
A. )C 


* W V> C> 

h 

lU liJ (J u 

to yj v> to 

3 


(L 

o e • 

* ♦ I 

z < 


- W - • N 


u o 

V S N. I 
N N J 

r Q u. 
a e r 
t ^ < 

C ‘u IL 
^ ^ u. 

^ 'Z 
<0 tn • I 

• • IL 


«-« 1^ 

X. 


u a 

N N 
U. 

X • 

<■ z 

2 S 

• n 


-• u • 

I N U ' 
•ON 


• ly N 

^ w* I 
• ^ ' 

X ^ 

111 tv 

< ^ 


u o 

V. • 

O M 


lU I 


Z ^ • V ^ 


* UJ 

• ^ 
X • 


a ■«* 


tA'V*. U ; 

3 3 O : 


- O ji :Z a. 


O 3 
:t to < 

O “ 


, in H 

z 

yj a 

dig 

W h U 


X a 
u o 

>-■ t- 
J < 

< z 


-J J 
a tL : 

n K3 ;5 -I 

p“ 1-» 

ir n. 

r a c. tr 1 

JJ ’3 ^ lU 
^ K K < 
J < < J 
« X z t 


o • , 

Z • _i J O I 

?v * n vj X a o 

<v _l tV ^ O o z 

« L. — 

3 ’i ct K ' 

a ■< (t :»■ a a < 

!y u ia a ‘J I 

a »- K ^ u I 

_i o J J < tf ly I 

•t < < < a X irt ' 


t p 

i a c. 
O 3 

z z 

' M a 

! <r < 

d a 


1 N 

o * 

p a 

• 3 

a * 

^ «y 

< V 


** a u 

S Z • N 

K >• ‘-I CV 

z r w 3 3 

N rj — s r 

^ < t i z Cl 

N - # • N 

X a J - 

N N *L -M 

< k. 


fc. • ?> > 


It IT 
< < 
ii a 


Z -J 
tc. a 


Z D 
3 

ly o c. 

> 4 
J >* o 

;i a a 

to X < 


li. ¥>'■ 
h' 

• Ul 

Z 3 ^ 

• a. 

> ^ V 
N N m s 

a w 3 • 

o 3 fK • 
Z Q t“ • 

• a • • 

2 • V ro 

• n o 3 
0 3 z a 
N a 3 a 
N X O 

a 2 

3T > ft: 

•IT 3 3 
< 3 J < 

^ uJ O X 


• z 

u • - 

N Ul N Z 
t* N tv ♦ 

• N Ul U 
Z X N 

U (V 


Q 3 I 
J K 
U- “ 

Z tv ' 

N 3 
• C 

H» N a 
3 

rj a 
w a _j 

? .• iL 

PIT 

a *y , 
a < n 

J a 3 
J < 0 
■a j- '*!- 
< n 3 
J X n 


\ V • N 


. u * 

! N ri 
PJ W 

■ a V 

> u • 
[ N CV 


I rv •* 

tsi 

k 

I . 2 

UJ a 


G < 
to X 


- 3 • 

L».- ^ CV 

u • a 

V a o 

• p • 

X a - «» 

tv X 
< 

« M » 2 ^ 

•X 3 o a * 

— a iL a ^ 

< a K 

[ii .1 

Z 7. Z ^ 

X a a _s pf _i 
a a a 'll ‘J *- < 

♦- I- K a K- « o 

< < 4 < J X z 

X X X J < u u 


Ul 


I, 

QUAxjtsj 


2-6 


Fig. 2-1 Rigid Format 7 ALTER-DMAP Statements (with Free Surface Displacement Recovery) 
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Fig. 2-2 Rigid Format 7 DIWAP Listing (with Free Surface Displacement Recovery) (Sheet 1 of 8) 
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Fig. 2-2 Rigid Format 7 DMAP Listing (with Free Surface Displacement Recovery (Sheet 6 of 8) 
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Fig. 2-2 Rigid Format 7 DMAP Listing (with Free Surface Displacement Recovery) (Sheet 8 of 8) 
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The following set of DMAP statements are utilized when gravitational stiffness 
is ignored and free surface displacement data is not desired. Copies of the ALTER- 
DMAP statements and the complete DMAP listing are presented in Figs. 2-3 and 2-4, 
respectively. 


136a. Partition 

area matrix 


[■^dd] 


to form the fluid inertance matrix j 

^m] 


[-da] = 


and the surface 


Ku 

0 

bp 

A 

0 

em 



136b. Calculate the transpose of , [a^^J 


136c. Solve for the pressure recovery matrix 

hus] - [%"] [^e®] 

136d. Calculate the fluid mass matrix 

[^fld] " [^em] fdu2] 
136e. Output into cards or tape 
170. EXIT. 


2. 2 DESCRIPTION OF DMAP OPERATIONS FOR HYDROELASTIC NORIVIAL MODES 
ANALYSIS 

The following is an outline of the major operations performed in Rigid Format 3 
for free- vibration analysm of a hydroelastie system. The special case of zero free 
surface pressure has been verified in test eases and the general case which has not 
been verified is also presented . The modifications required for a Rigid Format 13 
hydroelastic analysis with differential stiffness are nearly identical as those required 
in Rigid Format 3; thus a description of DMAP operations is not included. ALTER- 
DMAP statements and modified Rigid Format DMAP listings, however, are presented 
for Rigid Format 13 as well as Rigid Format 3. 
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Fig. 2-4 Rigid Format 7-DMAP Listing (without Free Surface Displacement Recovery) (Sheet 5 of 7) 
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3. GPl generates coordinate system transformation matrices, tables of grid 
point locations, and tables for relating internal and external grid point 
numbers. 


6. GP2 generates Element Connection Table with internal indices. 

8, PLTSET transforms user input into a form used to drive structure plotter. 
10. PRTMSG prints error messages associated with structure plotter. 

13i Go to DMAP No. 17 if no undeformed structure plot request. 

14. PIXDT generates all requested Undeformed structure plots. 

16 . PRTMSG prints plotter data and engineering data for each undeformed 


19. GP3 generates Grid Point Temperature Table. 

21. TAl generates element tables for use in matrix assembly and stress 
recovery. 

26. SMAl generates stiffness matrix [h^g] ^rid Point Singularity Table. 
28. SMA2 generates mass matrix ^^ggj • 

37. Equivalence [^^g] [^gg] general elements. 

39. Go to DMAP No. 42 if no general elements. 

40. SMA3 adds general elements to stiffness matrix [^^g] ^ Obtain stiffness 

mstrfct [Kgg] . 

44. GP4 generates flags defining members of various displacement sets 
(USET) and forms multipoint eonstraiht equations j |ug| ^ 0. 

46. Go to DMAP No. 118 and print error message if no independent degrees 
of freedom are defined. 

48. Equiv^ence j'KgJ to |:K^j and j'Mggj to |M^j if no multipoint 
constraints. 

50. Go to DMAP No. 54 if general elements present. 
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51. GPSP determines if possible grid point singularities remain. 

52. OFP formats table of possible grid point singularities and places it on the 
system output file for printing. 

55. Go to DMAP No. 60 if MCEl and MCE2 have already been executed for 
current set of multipoint constraints. 

56 . MG El partitions multipoint constraint equations [^^g] ” 

solves for multipoint constraint transformation matrix j*G^ ] ^ ~ [^m] ^ 

N- 

58. MCE 2 partitions stiffness and mass matrices 



k 

K 

nm 
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nn 

M: 

nm 


Vn 

K 

mm 

and 

Kg] ' 

mn 

M 

mm 


and performs matrix reductions 


[KJ - KV pm] ^ [°m] pmm] pm] 


and 


K] = [«»n] * [°m] [«mn] ^ [«L] M * pm] [“mm] pm] 

61. Equivalence an<^ [^nn] if no single-point 

constraints . 

63. Go to DMAP No. 66 if no single point constraints. 

64. SG El partitions out singler^point Constraints 


Kl' 

% 

K„ 

fs 

and [M^] - 


1 

fs 

1 

5f 
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ss 


^f 

L 

m/ 

ss 


67. Equivalence to jk^^j and jM^^j if no omitted coordinates 

69. Go to DMAP No. 74 if no omitted coordinates. 
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70. SMPl partitions constrained stiffness matrix 


N = 


K 


aa 


K 


ao 


K. 


oa 




oo 


solves for transformation matrix Tg 1 = - Fk 1”^ Tk 1 

L oj L ooJ L oaj 

T 

and performs matrix reduction ~ [^aa] ^ [^o] 

71. Form the total structure plus fluid gravitational and/or ullage stiffness 
matrix 

[^L] = Paa] * Kb] 

72. SMP2 partitions constrained mass matrix 


[%] 




aa 


M. 


ao 


M 


oa 


M 


00 


and performs matrix reduction 

Ka] - Paa] " Poa] [®o] " [°o] Poa] * [<] poo] Po] 

74. Form the total structure plus fluid mass matrix 

PU'PaaPPffld] 

75. Go to DMAP No. 84 if n® free-^body supports. 

76 . RBMGl partitions out f ree-body supports 


iv]" 



K; 

"/r' 





and pVI 1 = 




t aaj 


K 0 

K 



_ rjP 

rr 



M, 




M 


h 


M 


rt 

;78. RBMG2 decomposes constrained stiffness matrix 

P«] = [%] p«l 


M 


rr 
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80. RBMG3 forms rigid body transformation matrix 

[“] ' - [%] 

calculates rigid body check matrix 


M = Kr] 

and calculates rigid body error ratio 



82. RBMG4 forms rigid body mass matrix 

[”r] = [«rr] " [«fr] M * [“/r] " [“T [“w] [“]• 


85. . DPD extracts Eigenvalue Extraction Bata from Bynamics data block. 


87. Goto.BMAP No. 116 and print error message if no Eigenvalue Extraction 
Bata. 

89. REAB extracts real eigenvalues from the equation 

] iu 1 = 0 

[_ aa aaj ) a/ 

calculates rigid body modes by finding a square matrix [^j-o] 


such that 


KJ " [*ro] K] i'^ro] 


is diagonal and normalized, computes rigid body eigenvectors 


["^ao] 


ro 


<f> 


ro 


calculates modal mass matrix 

T 


H '['^alKa] M 
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and normalizes eigenvectors according to one of the following user 
requests; 

1) Unit value of selected coordinate 

2) Unit value of largest component 

3) Unit value of generalized mass. 

92. OFP formats eigenvalues and summary of eigenvalue extraction informa- 
tion and places them on the system output file for printing. 

94. Go to DMAP No. 120 and exit if no eigenvalues foimd. 

95. SDRl recovers dependent components of the eigenvectors 


(^f ] 

1*.| ■ W • (*.) ■ [”.] IM • 

fcl- W 


and recovers single point forces of constraint j * 

98. Equivalence SIL to SIP and BGPDT to BGPDP when one or more geometric 
grid points exists. 

101. PLTTflAN modifies BGPDT and SIL for functional modules SDR2 and 
PLOT. 


105. SDB2 calculates element forces and Stresses (OEFl, OESl) and prepares 
eigenvectors and single-point froces of constraint for output (OPHIG, 
PPHIG, OOGl). 

106. OFP formats tables prepared by SDR2 and places them on the system 
output file for printing. 
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107a. 


107b. 


107c. 


107d, 


10 7e. 


107e. 


107f 


107g. 


108. 

109. 

111 . 


Calculate the modal stiffness matrix 

K] = [*a] Ka] [*a] 

Solve 

[a,4] = [mj-1] [k^] 

Calculate the modal pressures (pressure deviations in the general case) 

[VJ = [^au.] K] r<] 

Calculate the free surface modal displacements 

[^l] “ [’’uia] [^a] 

Print . [^i] 

When recovery of free surface displacements is not desired (107d) and 
(107e) are replaced by the statement. 

Print ['<], [P/^] 

When the gravitational and/or ullage stiffness are included in the analysis 
the following additional operations are required for complete pressure 
recovery; 

Calculate the reduced out modal pressures 
Calculate the remaining modal pressures 

Go to DMAP No. 112 if no deformed structure plots are requested. 

PLOT genera.tes all requested deformed structure plots. 

PRTMSG priiits plotter data and engineering data for each deformed plot 
generated. 
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113. Go to DMAP No. 120 and make normal exit. 

115. NORMAL MODE ANALYSIS ERROR MESSAGE NO. 1 - MASS MATRIX 
REQUIRED FOR REAL EIGENVALUE ANALYSIS. 

117. NORMAL MODE ANALYSE ERROR MESSAGE NO. 2 - EIGENVALUE 
EXTRACTION DATA REQUIRED FOR REAL EIGENVALUE ANALYSIS. 

119. NORMAL MODE ANALYSIS ERROR MESSAGE NO. 3 - NO INDEPENDENT 
DEGREES OF FREEDOM HAVE BEEN DEFINED. 


Copies of the ALTER-DMAP statements and the complete DMAP listing for the 
Rigid Format 3 hydroelastie normal modes analysis are presented in Figs. 2-5 and 
2-6, respectively. The corresponding information for the Rigid Format 13 hydro- 
elastic normal modes analysis with differentia’ stiffness is presented in Figs 2-7 and 
2-8, respectively. 

2.3 PROBLEM DEFINITION AND BULK DATA CONSIDERATIONS 


2. 3. 1 Grid Sequence 


The modified NASTRAN hydroelastie technique defined by the altered Rigid 
Formats requires special purpose bulk data for individual problem definition. The 
Rigid Format 7 grid set must be ordered according to the following convention: 

(1) generalized free sirnface pressures, |p^| 

(2) generalized structural interface pressures, 

(3) generalized internal pressures, |pj| 

' (4) discrete free surface displacements, 

(5) harmonic free surface displacements, 

(6) discrete structural displacements, 

(7) harmonic structural displacements, |U„ ) 

I ®h / 


{"4 
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Fig. 2^ Rigid Format 3 DMAP Listing (Sheet 4 of 6) 
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Volume n 


When free surface displacement recovery is not desired, categories (4) and (5) 
are absent. Although the external grid indexing* Ig, is such that the generalized 
pressures (1), (2) and (3) are of lowest indices, the NASTRAN internal indices for 
the pressures Ij, are; 

Ij = 1, 000, 000 X + Ig 

where 

I = N+i for a cosine series 
n 

I = N+1/2 for a sine series 
n 

In order to avoid conflicting indices and to accommodate partitioning operations, 
ejfternal grid index bounds must be strictly adhered to (Note; external and internal 
indices for displacements are equivalent). The chosen external index bounds are: 

1 llg (P 3 , Pf. Pi) < 999 

lO&O < 4 (U,, U. ) < 2999 
E f % 

3000 < Ig (Ug) < 9999 


10*000 < I„ (Uq ) < 999, 000 
E 


The above sequencing bounds also hold for Ug 
for modal hydroelastic analysis. 

2. 3. 2 Coordinate Systems 

The most convenient displacement set coordinate systems for hydroelastic 
problems are local spherical or cyliindrical reference frames in Which one of the linear 
displacements is norm^ to the fluid bounding surface and the other two linear displace- 
ments are tangent to bounding surface; such reference frames are also ideal for 
structural grid points defining the surface of a plate or shell not interfacing with a 
fluid. In the case of the free surface. When displacement recovery is desired, a 
cyMndrical reference frahie with the ig (z) displacement eomponenfe directed normal 
to the free surface should be defined. 


and Ujj in Rigid Formats 3 and 13 
®h 
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Typical referettse frames with defining parameters are illustrated in Fig. 2-9. 
2. 3. 3 Constraints and Reductions 


The required and recommended constraints in the EUgid Format 7 fluid data 
calculation and Rigid Format 3 or 13 hydroelastic normal modes analysis follow 
differing rules due to the special purpose of each respective calculatiGn. The final 
analysis sets, however, must be compatible. The definition of local spherical and 
cylindrical reference frames described in Section 5. 3. 2 plays a crucial role in the 
considerations outlined below. 

The constraints and reductions applied in the Rigid Format 7 fluid data caleula- 
tipn are: 

• Multipoint Constraints (MFC) - The multipoint constraints are used for 
harmonic transformation of discrete free surface and structural displace- 
ments. Only the displacement components to be retained in the analysis set 
should be transformed (e.g. , outward normal only, normal and tangential 
only). As a general rule only those free surface displacements normal to 
the surface are transformed. 

• Single Point Gonstralnts (SFC) - The single point constraint set includes all 
structural displacements (Harmonic or discrete) to be omitted by Guyan 
reduction in the Rigid Format 3 or 13 modal analysis; components normal to 
the fluld/structure interface should not be constrained. It should be noted 
that there is no structural stiffness matrix to be decomposed in the Rigid 
Format 7 calculation. Single point constraints defining physical structural 
constraints are included in the SPC set. 

In the cases where fluid free surface displacement sets are not defined, the 
free surface pressures becom.e members of the SPG set; in addition, for 
intermediate fluid fill conditions the pressures above the intermediate free 
surface become members of the SPC set. 

• Guyan Reduction (OMIT) The OMIT set must include all fluid Internal 
pressures. Strueturai displacements may not appear in this set since the 
structural stiffness matrix is not defined in Rigid Format 7. 
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• Special Partitioning -• A partitioning vector, is defined by DMI state- 

ments to partition the analysis set into generalized pressure and generalized 
displacement subsets. If fluid free surface displacements have not been 
defined in the formulation, further partitioning is not required. 

When fluid free surface displacements are defined partitioning of the pressure 
and displacement subsets is required. The vector |g 9| is defined by DMI 
statements to partition the pressure subset into free surface pressure and 
strucmral interface pressure subsets. The vector is defined to parti- 

tion the displacement subset into free surface and structural displacement 
subsets. It should be noted that the structural displacement subset corre- 
sponds to the Rigid Format 3 analysis set to be defined below. 

The constraints and reductions applied in the Rigid Format 3 or 13 hydroelastic 
normal modes analyses are: 



I ' I 

L 1 ! 



i i : 
L J i 


I } i 

( i • 

t 

i 



• Multipoint Constraints (MPC) ^ The multipoint constraints are used for 
harmonic transformation Of discrete structural displacements if harmonic 
reduction is utilized . All six displacement degrees of freedom must be 
transformed for each grid point associated with the harmonic constraints. 
Additional multipoint constraints not associated with the harmonic transfor- 
mation may also be defined. 

• Single Point Gpnstraints (SPC) Single point constraints are utilized to define 
structural supports on harmonic as well as discrete displacements as 
required . 

• Guy an Reduction (OMIT) r - The OMIT set includes harmonic displacement 
components to be reduced out (e.g, , tangential displacements and rotations). 
Discrete displacement components are also members of this set as required. 

2. 3. 4 Special Bulk Data Considerations 

The bulk data required in the modified Rigid Format 7 calculation of fluid matrix 
data closely follows the unmodified Rigid Format 7 bulk data with a few exceptions. 
These are; 



r I j 




(1) A free surface may not be defined in terms of FSLIST. If free surface 
displacement recovery is required (or gravitational stiffness is included) 






Volume n 


the free surface is defined as part of the BDYLIST, If free surface 
displacement recovery is not desired the BDYLIST only includes the fluid 
points at the structural interface. 

(2) A dummy EIGC card is required to avoid an error message. 

(3) A PARAM card is required with the following data: 

Field 1 = PARAM 
Field 2 = A999 

Field 3 = N (floating point number) 

Field 4=0. 

Where N is equal to the number of symmetric sections of structural boundary 
about the fluid circumference being modeled by structural elements. 

(4) The partitioning vectors t |g 9| and |up| must be defined on DMI 

cards as described in Section 5.3. 3. 

The fluid matrix output from the Rigid Format 7 run may be placed on cards o:’ 
tape. If it is output onto DMI cards, these cards must be Used as part of the input 
data for a Rigid Format 3 or 13 hydroelastic normal modes analysis run. The 
remaining bulk data for the hydroelastic normal modes analysis is the same as for a 
dry structure. 
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3 - 1/8-SCALE SPACE SHUTTLE EXTERNAL TANK HYDROELASTIC ANALYSIS 


An investigation of 1/8-scale external tank dynamics in a free-free supported 
condition utilizing the modified versions of Rigid Formats 7 and 3 has been under- 
taken. The finite element hydroelastic model illustrated in Figs. 3-1 through 3-4 is 
described by a grid set consisting of 348 pressure degrees of freedom, 2,058 struc- 
tural degrees of freedom, and 768 harmonic structural degrees of freedom. The 
structural model contains elements representing all rings, stringers, internal 
stiffening members and asymmetries present in the 1/8-scale test article. Normal 
and tangential motion of the tank surface is conveniently described in terms of local 
spherical and cylindrical reference frames. 

The fluid model deseribing the dynamics of the fluids consists of 204 and 144 
degrees of freedom for the LOX and LH 2 model fluids, respectively, in the liftoff 
condition; harmonics n == 0, 1, 2, 3 were chosen to describe asjmmetric dynamics 
with the pitch plane taken as an axis of symmetry. The respective LOg and LHg free 
surface motions in this application were not desired. Free surface displacement sets 
are therefore not defined, however the fluid free surface condition is described in 
Rigid Format 7 by application of single point constraints (SPC's) on free surface 
pressures. In addition, for liquid levels other than the liftoff condition, single point 
constraints are applied to pressure degrees of freedom at and above the appropriate 
free surfaces. This feature of the analysis provides capability for study of many 
liquid fill conditions utilizing a single grid set fluid model. 

Harmonic reduction retaining harmonies n = 0, 1, 2, 3 was applied to the 
structural grid of the external tank (ET) for consistency with the fluid pressure grid. 
The multipoint constraint data for this study was generated by a general purpose data 
preparation program HARM which is discussed in Appendix A. Since Consistency is 
desirable mainly at the fluid /structure interface, the harmonic transformation was 
limited to the interface points associated with the liftoff condition. Many of the 
remaining grid set displacements were ’’omitted" by Guyan reduction. In addition, in 
the basic analysis, all generalized rotations and tangential displacements in the 
harmonic set of displacements were "omitted". The analysis set of displacements 
therefore consists of 128 harmonic, outward normal displacement degrees of freedom 
(32 meridional locations, harmonics n = 0, 1, 2, 3). 
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Three liquid fill conditions have been studied consisting of liftoff, post max Q 
and empty. In terms of liquid height above the respective bulkheads the conditions 
are identified as: 






liftoff: 

post max - Q: 


= 190.5 cm, hTu = 
LH2 

^LOX ~ ^ 


= 358 . 14 cm 

, „ = 330.2 cm 
L.M2 


. empty: = 0 

The liquid free surfaces are normal to the tank axis since the NASTRAN hydro- 
elastic analysis is currently limited to axisjrnimetric fluid configurations. 


Normal modes analysis for each of the fill conditions was performed with the 
modified version of Rigid Format 3. All 128 eigenvalues and 25 mode shapes and 
modal pressure distributions were extracted by the Given's method. Computational 
efficiency was quite satisfactory for this rather complex mathematical model. Com- 
putation time on the Grumman IBM 370/165 computer for the fluid mass and pressure 
recovery matrices (LOX and LK 2 combined) on Rigid Format 7 was on the average 
slightly less than 5 CPU minutes per liquid level and computation time for normal 
modes analysis in Rigid Format 3 was on the average about 18 CPU mmutes per liquid 
level (16 min. computation, 2 min. plot preparation). In previous attempts to study 
the d 3 Uiamics of the same finite element representation with the standard NASTRAN 
hydroelastie analysis, computation times were in excess of 70 CPU minutes with only 
one eigenvalue and mode shape extracted by the unsymmetrie inverse power method. 
Actually the computation time for normal modes analysis can be significantly reduced 
by utilization of NASTRAN's checkpoint- re start capability. By utilization of check- 
pointing in a first case liquid level NASTRAN run, structural and fluid data may be 
saved on tapes eliminating much of the computation time associated with preparation 
of matrix data common to all liquid levels (especially reduction of the structural model 
representation), total computation time per liquid level may be reduced to about 10 
CPU minutes for cases subsequent to the initial run. 

Modal data for each of the three fill conditions studied are summarized in 
Tables 3-=-l throu^ 3-3 with dome pressure gain presented as a measure of relative 
significance in Pogo susceptibility. The pressure gain is defined as the dome pressure 
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(LO 2 and LH 2 ) in the vicinity of the feedline interface resulting from a unit modal 
acceleration. Plots of the mode shapes are presented in Figs. 3-5 throu:rh 3-8. 
Unrealistic behavior of the individual tank dome apexes is present in the liftoff con- 
figuration illustrated in Figs. 3-5 through 3-26 and in the post max Q configuration 
illustrated in Figs. 3-27 through 3-31. This localized behavior was eliminated by 
incorporation of an "apex fix" in the post max Q configuration illustrated in Figs. 3-32 
through 3-53 and in the empty configuration illustrated in Figs. 3-54 through 3-75 
with negligible effect on natural frequencies and overall mode shapes. The "apex fix" 
consists of a set of multipoint constraints (MFC’s) forcing the apex and the set of 
grid points connected to it by triangular elements to move as a rigid body in the 
zeroth harmonic as shown in Fig. 3-76. 


An assessment of the accuracy of the current NASTRAN analysis results was 
made on the basis of available test and previous analysis data. Comparisons of the 
current results with limited test, previous NASTRAN finite element, and NASTRAN 
equivalent beam model results are presented in Table 3-4. For the liftoff condition, 
the current and previous NASTRAN first axial mode frequencies are in disagreement. 
This is attributed to the extreme Guyan reduction required in the previous analysis to 
keep computation time at a reasonable level. Test data for the liftcff condition is not 
available at the present time. In the case of the post max Q fill condition, limited 
test results are available with some uncertainty in the second overall bending fre- 
quency. Excellent agreement between current NASTRAN finite element and test 
frequencies in the first axial mode and poor agreement in the reapeetive bending 
frequencies were noted; the bending mode discrepancy is also present in the empty 
condition results. An effort to locate the source of discrepancy was initiated. The 
present first bending mode and n = 3 shell mode of the cylinder presented in 
Figs. 3-54 and 3-55 respectively, compared well with the previous NASTRAN results 
(see Table 3-4). No comparisons were made with the third empty External Tank mode 
calculated in the previous analysis since this mode is predominently an n = 4 shell 
mode and the current analysis is limited to the harmonics n = 0, 1, 2, 3. 

A comparison of the current NASTRAN and test firing bending frequencies in the 
empty condition indicates a frequency error 


U) Test Empty ^ 89 
to Anal Empty 105.6 


(3-1) 
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which applied as a correction to the post max Q first bending frequency results in 

It) corrected = 0. 8428 x 49. 2 Hz = 41. 5 Hz (3-2) 

which Is much closer to the test result, 37 Hz. This leads one to suspect strongly 
that an inaccuracy exists in the structural finite element model. A series of modifi- 
cations in the structural model were Implemented in search of sources of error. In 
one case, the large elements on the LH 2 tank cylinder were converted from CQUAD2 
membrane -pi ate elements to CQDMEM membrane elements since it was suspected 
that individual plate bending strain energy was unrealistically high; the results of the 
study showed no change in the fundamental bending mode frequency. 

In another attempt to find the source of inaccuracy, ring frame torsional and 
axial stiffnesses were reduced to reflect partial effectiveness in the axial direetlon; 
again there was no change in the fundamental bending mode frequency. At this point 
the effort to locate the discrepancy was halted because further effort was considered 
to be outside the scope of the present task. 

In order to verify the fluid model, modal pressure distributions in the first 
axial and bending modes for the post max Q fill conditions were studied; the pressure 
distributions are illustrated in Figs. 3-77 and 3-78. Since nearly all of the system 
mass is associated with the fluids, the pressure distributions in the free-free 
supported modes must reflect an overall state of equilibrium. This condition is 
apparent in the axial (n = 0) pressure distributions by the nearly equal and opposite 
bulkhead pressures {note that a small n = 0 pressure eontribution appears in the 
bending mode since the tank structure is not purely axisymmetric). In the bending 
mode the net lateral force and pitch moment are nearly zero; thus confidence is 
gained in the fluid representation in the lower mode shapes. Further verification of 
the fluid representation in the higher modes is effected by examination of pressure 
distributions in the LOX tank in the post max Q fill condition illustrated in Fig. 3-79; 
it is noted that the pressure and displacement profiles below the level of the free 
surface have similar shapes with nearly coincident vibration nodes. In addition, 
pressure distributions near the free sirnface in all lateral pressure distributions 
(n > 1) approach zero pressure indicating the "effective slosh mass" contributions 
which are out of phase with shell motion in the hydroelastic frequency range. Further 
quantitative checks were made by hand computation of modal generalized mass of the 



Volume n 


fluid on the basis of pressure distributions which compared well with approximated 
modal mass values based on a NASA LaRC beam model and with the NASTRAN 
computed modal mass. With verification of the fluid representation fairly complete, 
a checkout of the structural model was initiated. 

An initial verification of the aceuracy of the harmonic reduction scheme applied 
the structurally asymmetric tank was made by comparison of the current empty tank 
results with previous NASTRAN results for the same finite element model without any 
reductions. Only three modes were obtained by the Inverse power method in the 
previous analysis and these are presented in Figs. 3-80 through 3-82. 

In the event that the reader wishes to review and/or revise the l/8-seale 
external tank mathematical model, he must be aware of the fact that some shell grid 
point displacements are not expressed in terms of local cylindrical or spherical 
reference frames. Thus displacements associated with shell grid points 8134, 8352, 
8355 must not be related to the harmonic displacement set; it should also be noted 
that any grid point displacements not associated with shell surface motion may not be 
related to the harmonic set. 
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liable 3-2 1 / 8 -Scale {External TankiHydroelastic Mode Summary {at Liftoff) 
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Table 3-3 Empty 1 /8-Scale 


Mode No. 

Freq. IMz) 

Modal II 

4 

105.6 

0.0425 

5 

153.0 

0.0094 

6 

161.7 

0.0172 

7 

226.0 

0.0497 

8 

257.8 

0.0770 

9 

274.7 

0.0271 

10 

328.3 

0.0122 

11 

332.0 

0.0149 

12 

332.8 

0.0234 

13 

343.7 

0.0118 

14 

357.8 

0.0696 

15 

431.0 

0.0210 

16 

459.1 

0.0615 

17 

472.9 

0.0114 

18 

482.2 

0.0185 

19 

498.6 

0.0076 

20 

513.2 

0.0697 

21 

533.1 

0.0243 

22 

567.2 

0,0487 

23 

604.6 

0.0391 

24 

625,4 

0.0116 

25 

628.0 

0.0144 


Tank Mode Summary 




Description of Mode 


ET 1st Bending n=1 
LH 2 Cylinder n=3 
LH 2 Cylinder hr2,3 
ET 2nd Bending n=1 
ET 1st Axial n=0 
LM 2 Cylinder n=2, 3 
LH2 Cylinder, LOX n=3, 2 
LOX. LH2 n=3,2 
LOX, LH2 n=3,2 
ETn=3 

ET Bending n=1,3 

LM2 Cylinder n=2 

LM 2 Cylinder n=3, ET n=1 

LM 2 Cylinder n-3 

ETn=2 

LOX n=3 

ET n=1 , 2, 3 

Et n=2, 1 

ETn=1,2,3 

ET n=2, 1 , 3 

LOX n=3 

Et n=3,2 









Volume II 



3-14 




f 









Volume II 



3-15 


Fig. 3-7 Liftoff Mode 6 
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Fig, 3-8 Liftoff Mode 7 
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Fig: 3-9 Liftoffi Mode 8 
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Fig. 3-10 Liftoff Mode 9 
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Fig. 3-11 Liftoff Mode 10 
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Fig. 3-12 Liftoff Mode 11 
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Fig. 3-13 Liftoff Mode 12 
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Fig. 3-14 Liftoff Mode 13 
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Fig. 3r16 Liftoff Mode 15 
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Figi3-17 Liftoff Mode 16 
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Fig; 3-18 Liftoff Mode 17 
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Fig^3-19 Liftoff Mode 18 
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Fig; 3-20 Liftoff Mode 19 
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Fig. 3-21 Liftoff Mode 20 
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Fig. 3-22 Liftoff Mode 21 
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Figi3-73 Liftoff Mode 22 
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Rg. 3-24 Liftoff Mode 23 


Volume n 



I 

3-33 I 


Fig. 3-25 Liftoff Mode 24 
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Fig. 3-26 Liftoff Mode 25 



Volume II 



3-36 


Fig, 3-28 Post Max Q Mode 5 {without Apex Fix) 



Fig: 3-29 Post Max Q Mode 6 (without Apex Fix) 
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Fig. 3-31 Post Max Q Mode 15 (without Apex Fix) 




Fig; 3-33 Post Max Q Mode 5 
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Fig. 3-36 Post Max Q Mode 8 
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Fig. 3-38 Post Max Q Mode 10 
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Fig. 3-39 Post Max Q Mode 11 
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Fig. 3-40 Post Max Q Mode 12 



Volume n 



3-50 




Volume II 



3-53 


Fig. 345 Post Max Q Model? 
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Fig. 3^3 Post Max Q Mode 25 



Volume II 



3-63 




Fig. 3*55 Empty External Tank, Mode 5 



Fig: 3-56 Empty External Tank, Mode 6 
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Fig. 3-58 Empty External Tank. Mode 8 
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Fig. 3^59 Empty External Tank, Mode 9 
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Fig. 3-61 Empty External Tank, Mode 11 
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Fig; 3-62 Empty External Tank, Mode 12 




Volume n 



Volume n 



3-73 


Fig. 3-65 Empty External Tank, Mode 15 
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Fig4 3^66 Empty External Tank, Mode 16 
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Fig. S-S? Empty External Tank, Mode ‘•7 
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Fig. 3^8 Empty External Tank, Mode 18 
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Fig. 3^9 Empty External Tank, Mode 19 
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Fig^ 3-71 Empty External Tarik, Mode it 
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Fig^ 3*73 Empty External Tank, Mode 23 



Fig; 3-74 Empty External Tank, Mode 24 
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Fig. 3*79 Two Post Max Q External Tank Modes with Modal Pressure (LOo Tank Illustrated) 





Fig, 3^0 Em^ity External Tank, Mods 1 list Bending) Calculated with No Reductions 
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APPENDIX A - ’’HARM" - A FORTRAN PROGRAM FOR AUTOMATIC 
HARMONIC MULTIPOINT CONSTRAINT DATA GENERATION 


The FORTRAN IV computer progrmn listing presented in Fig. A-i has been 
written to aid the analyst In preparation of harmonic transformation data required in 
Rigid Format 7 fluid matrix data csdculations and in Rigid Format 3 or 13 hydro- 
elastic normal modes calculations. The harmonic transformations are limited to 
descriptions of symmetric djTiamics of geometrically ajdsjrmmetric bodies. The 
author has utilized this program on the Grumman Calldata time sharing system. On 
this system h^ut data was prepared on a file, FTO4F001 and the ou^ut data was 
placed on a file FTO7F001 {note the read and write parameters 4 and 7 respectively, 
in the listing). The output file is punched onto cards and listed in typical runs. 

The input data for "HARM" is as follows : 

NH = number of harmonics 

NPT = number of grid points to be transformed 

NFS = parameter ^ 

if = 0 transform only the "1" degree of freedom (the radial) 

and 

if 0 transform all 6 degrees of freedom. 

KHARM ^), 1 = 1, NH; harmonic numbers to be included. 

For each of the "nPT" grid points, the following li^ut data is required: 

NGR\ = grid point identification number 

NROW = meridianal row number associated with the grid point; convention 
chosen by analyst 

THET = circumferential location of grid point in degrees 

NTYP = parameter identifying reference frame type 
= 1 (cylindrical reference frame) 

= 2 (local spherical reference frame) 

It should be noted that harmonic grid point identification numbers, KH (I), 
following the convention KH (I) = 1000000 + 100 * NROW + KHARM (I) are generated by 
"HARM", The analyst Should separately prepare GRID or GRIDB cards for NASTRAN 
bulk data. 
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Fig. A-1 “HARM'' Fortran Listing (Sheet T of 4| 
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